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Abstract: The present study involves the CFD analysis for the 
prediction of swirl effect on the characteristics of a steady, 
incompressible flow through an S-shaped diffusing duct The 
curved diffuser considered in the present case has S-shaped 
diffusing duct having an area ratio of 1.9, length of 300 mm 
and turning angle of 22. 5 c V22.5°. The static pressure, total 
pressure, velocity and turbulence intensity were accounted. 
The improvement is observed for both, clockwise and anti- 
clockwise swirl, the improvement being higher for clockwise 
swirl. Flow uniformity at the exit is more uniform for 
clockwise swirl at the inlet. 

Keywords: Curved diff users, intake ducts, swirling flow, 
secondary flows, pressure recovery 

INTRODUCTION 

Diffusers are one of the standard challenges in fluid mechanics. 
The task of a diffuser is to decelerate the flow and to regain total 
pressure. It is more difficult to arrange for an efficient 
deceleration of flow than it is to obtain an efficient acceleration 5 . 
There is a natural tendency in a diffusing process for the flow to 
break away from the walls of the diverging passage reverse its 
direction and flow back in direction of the pressure gradient. If 
the divergence is too rapid, this may result in the formation of 
eddies with consequent transfer of some kinetic energy into 
internal energy and a reduction in useful pressure rise. 

LITERATURE REVIEW 

The flow development in curved ducts with inflexion in 
curvature without diffusion has been reported by Rowe 1 
[1970]. The experiments were carried out in S-shaped ducts of 
mild curvature with the angle of turn being 45745° and 30730°. 
He showed that in a 45745° S-bend in which the secondary flow 
causes a complete interchange of wall fluid and fluid in the 
central core, possibly indicating the reason why heat transfer 
effects are large in wavy pipes. 

Bradshaw 2 [1972] have studied the flow characteristics within 
an S-shaped duct with strong curvature. They measured only 
pressures and observed a pair of counter-rotating vortices 
at the exit. The ideal pressure recovery coefficient is derived 
directly from basic principles. 

A computational investigation of flows in diffusing S-shaped 
ducts was provided by R.Menzies 3 . A high and low mass 
flow cases were compared with the experiments and 
calculations. The low mass flow case showed good 



comparison with experiment and calculations. A problem of 
the low mass flow case is the prediction of secondary flow. It 
was found that the SST turbulence model best predicted this 
feature. 

Anand 4 [2001] studied on the effect of the turning angle on 
the flow and performance characteristics of long S-shaped 
circular diffusers. He established the effect of the turning 
angle on the flow and performance characteristics of long S- 
shaped circular diffusers (length-inlet diameter ratio, L/Di = 
11.4 having an area ratio of 1.9 and centre-line length of 600 
mm.The experiments are carried out for three S-shaped 
circular diffusers having angles of turn of 15715°, 22.5722.5° 
and30730°. 

A detailed experimental study was been carried out by 
Whitelaw and Yu 5 [1993] on a scaled model of the RAE 2129 
S-shaped intake-diffusing duct having circular cross-section. 
The velocity and pressure characteristics of an S-shaped 
diffusing duct was measured in a water tunnel with three 
asymmetric inlet conditions using laser-Doppler anemometry 
at a Reynolds's number of 40,000. The three inlet conditions 
were introduced by thickening the boundary layer at various 
locations along the inlet circumference and corresponding to a 
side-slip, zero angle of attack and low angle of attack. S- 
shaped duct is considered as an intake passage for jet aircraft. 
They found out that a presence of enlarged vortices had 
impaired the non-uniformity of the corresponding stream wise 
mean velocity distribution at the exit. The maximum 
secondary flow velocity however, did not exceed 15% of the 
inlet bulk velocity in all cases investigated. 

CFD ANALYSIS 

Computational Fluid Dynamics (CFD) can be described as the 
use of computers to produce information about the ways in 
which fluids flow in given situations. CFD embraces a variety of 
technologies including mathematics, computer science, 
engineering and physics, and these disciplines have to be 
brought together to provide the means of modeling fluid flows. 
Such modeling is used in many fields of science and engineering 
but, if it is to be useful, the results that it yields must be a 
realistic simulation of a fluid in motion. 

PROPOSED WORK 

Ansys fluent is widely used for CFD analysis of various 
geometries. But In this study S-shaped diffuser is model 
created in ANSYS Design Modeller with thestructured multi- 
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block mesh generation in ICEM CFD. 

Simulations will be carried out in ANSYS FLUENT and 
appropriated boundary conditions for different faces were 
declared. The present study involves the CFD analysis for the 
prediction of swirl effect on the characteristics of a steady, 
incompressible flow through an S -shaped diffusing duct. The 
static pressure, total pressure, velocity and turbulence intensity 
were accounted. 




Fig. 1: Solid model of S -shaped diffusing duct 
Table 1: Dimension of the model 



Inlet diameter 


5 2.5mm 


Outlet diameter 


72.5mm 


Area ratio 


19 


Length of centre line 


300mm 


Offset 


5S.2mm 


Tmiiing angle 


22.5722.5° 


Radius 


3 S 1 nun 



The model is meshed in ICEM CFD. 



o Fluid was specified as air for the continuum type and the 

mesh was exported to Fluent for post processing. 
(STEP 2) Post Processing (In ANSYS Fluent): 

o Grid was checked and scaled. 

o 3D solver and segregated solution method was 
chosen. 

o Material property: Air was chosen as the fluid for 
flow, and its properties like density and viscosity at 
25 °C were selected. 

o k-co SST model is selected for turbulence with Low 
Re-corrections with the first cell y+ being maintained 
less than or equal to 1 at all wall locations . 

o At air inlet section, the inlet velocity of 27 m/s with 
different swirl intensity was specified. 

o Turbulence intensity of 3% based on inlet flow 
diameter was specified. At the exit section, the 
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Numerical Modeling Approach 
Simulation Procedure 

(STEP 1) Modeling (In ANSYS Design Modeller): 

Diffuser geometry is created by ANSYS Design Modeller using 
below dimensions 




Fig. 2: Meshed model of test diffuser 

o The model has been meshed with quadratic-mesh. 

Fine meshing with spacing 0.07 was done, mesh 

node's 245964 and mesh elements 229824. 
o Boundary conditions taken were for velocity at inlet, 

pressure at outlet and wall type, 
o Three simulations will be performed as per the 

experimental analysis with velocity magnitude of 

27m/s: 

1 . Uniform flow at the inlet 

2. Clockwise flow at the inlet (axial component = 
cos30°; tangential component^ sin 30°; radial 
component=0) 

3. Anti-clockwise flow at the inlet (axial 
component = cos30°; tangential component^ -sin 
30°; radial component = 0) 

The operating conditions will be kept at atmospheric 
pressure with the outlet being exposed to it, i.e., pressure 
outlet. All the walls will be considered as smooth no-slip 
walls. 

pressure was specified being equal to atmospheric 
pressure. 

o Second order upwind scheme was selected to solve 

continuity and momentum equations, 
o Convergence criteria of 10" 6 were taken, 
o Solution was initialized at inlet and made to iterate 

until it converges. 

Once solution is converged, various data for pressure and 
velocity were obtained and graphs were plotted. 

RESULT AND DISCUSSION 

Investigation of S -shaped diffuser has been carried 
out with help of FLUENT, a CFD tool to simulate the effect 
of swirl on the pressure recovery. The area ratio is kept 
constant at 1.9. The velocity profile was of 27m/s with 
different swirling intensity as 0°, 9°, 18°, 
22.574.5°,22.5713.5° and 22.5722.5°. 
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Clockwise flow chart 



Fig 3. Model showing planes and sections through which 
the parameters will be measured 

Analysis gives the effect of swirl angle at flow 
separation and on pressure recovery coefficient. 
Graphs of Flow Separation 

1) Variation of Static Pressure Recovery Coefficient and 
Total Pressure Loss 

Table 3: Showing the values of C p and C p .i oss for uniform 
and swirled flows at different sections 





uniform flo^v 


Clockwise-flow- 
30-deg 


Anti-clockwi&e- 
rlow-30-deg 


Posirion- 
Section 


Cp- 

uniform 

How 


uniform 
llow 


Cp- 

CW- 

30-deg 


Cp4055- 

nv-30- 


Cp-A- 

CW- 

30-des 


Cp4055- 

A-CW- 
30-des 


Section-0 


0.0S6 


0.0185 


0.0S6" 


0.0221 


0.0S6" 


0.022 1 


Section-9 


0.232 


D.D4TS 


0.2234 


D.0~09 


02230 


0.0~09 


Section- 
IS 


0.325 


0.0791 


0.3252 


0.1078 


0.3252 


0.1073 


Section- 

223/4.5 


0.:99 


0.1610 


:-a:~: 


0.1315 


2-.4I"2 


0.1315 


Section- 

22.5/13.5 


D.3S9 


0.1729 


0.4679 


0.1547 


0.46~9 


0.1543 


Section- 

22 5 '22 5 


3.487 


0.16S3 


0.5123 


0.1TSS 


0.5123 


0.1TSS 



Anti-clockwise-flow chart 




o 
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Fig 4. Variation of performance parameters along 
diffusing duct 



It is clear from the graph of Fig. 4 for a uniform flow at inlet: 
A) Uniform Flow Characteristics 



1. 



Static pressure recovery coefficient: In the first 
bend, the mass-averaged static pressure recovery 
coefficient increases continuously and then 
suddenly gets reduced slightly near to the point 
of inflexion due to change of curvature. The 
recovery coefficient again picks up just after the 
inflexion due to flow diffusion as a result of 
increase in cross-sectional area. 
Total pressure loss coefficient: The total 
pressure loss coefficient follows the same trend 
as that of the static pressure recovery coefficient. 
The overall static pressure recovery and total 
pressure loss coefficients are 0.487 and 0.1683 



of the inlet dynamic pressure 
respectively. 



(0.5pU z long ) 



The ideal pressure recovery (1 - 1/AR ) for this 
diffuser should be 72.3%. The lower pressure recovery 
compared with ideal static pressure recovery achieved may be 
due to the loss of inlet energy to creation of secondary 
flows and to overcome the higher surface resistance. 
B) Swirl Flow Characteristics: 
It is clear from the graph of Fig.4 for a swirl flow at inlet: 

1. Static pressure recovery coefficient: The mass- 
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averaged static pressure recovery coefficient 
increases gradually without getting affected by the 
inflexion in curvature. It is also observed that static 
pressure recovery coefficient for all the cases are 
comparatively much lower than the ideal recovery 
coefficient for this diffusing duct. 

2. Total pressure loss coefficient: The overall total 
pressure loss coefficient also increases for the swirl 
flow, which indicates high shear. 

3. The overall static pressure recovery and total 
pressure loss coefficients for both clockwise and 
anti-clockwise swirl are 0.5123and 0.1755 of the 
inlet dynamic pressure (0.5pU 2 i ong ) respectively. 

The increase of static pressure recovery could 
be due to the suppression of self -generated secondary 
flows as a result of swirl flow. It is observed that for both 
clockwise and anti-clockwise swirl flow, the overall static 
pressure recovery coefficient (-0.5123) is higher than that 
for uniform flow (-0.487). The slightly higher recovery 
for clockwise/anti-clockwise indicates that the loss in 
energy for the uniform flow case is slightly higher as a 
result of swirling angles. Thus, in all the above cases the 
low value static pressure recovery coefficient as 
compared to the ideal pressure recovery coefficient 
(-0.723) could be due to the combined effect of curvature 
resulting in the generation of secondary flows and surface 
roughness of the diffusing duct walls. 

Longitudinal velocity distribution for uniform flow 

Figure 5a gives the normalized longitudinal velocity 
(£4)ng/^ave(in)) distribution at different sections for uniform 
flow. Figure 4a (i) confirms that the fluid flow at inlet section 
is almost uniform except near the wall. At section 2, the fluid 
core shifts towards the outer surface under the action of 
centrifugal force. This initial distortion is observed up to 
section 4 (Fig. 5a (ii)), i.e., just downstream of the point of 
inflexion. There is no abrupt change observed between 
sections 2 and 3, and hence figures for those sections (section 
2 & 3) are not presented here. The bulk flow retains its 
position in the second bend. However, there is a change in the 
direction of curvature: the maximum velocities are observed 
near to the convex surface and minimum near the concave 
surface of second bend. No further qualitative changes are 
observed in the velocity distribution until the exit (Fig. 5a (iii)) 
of the diffusing duct. At the end of the first bend, stream-wise 
separation occurs very close to the inner surface of the first 
bend. The separated flow is seen from section 4 (Fig. 5a (ii)) 
and extends nearly up to the duct exit. 

Cross-flow velocity distribution for uniform flow 

Normalized cross-flow velocity ( C/ sec /C/ ave (in)) vectors 
for different planes of the duct are given in Fig. 4b. It is seen 
that even at the inlet (Fig. 4b (i)) very weak secondary flows 
are observed which could be seen as an upstream effect of the 
downstream curvature. At downstream stations this intensity 
increases and a pair of secondary vortices is observed dividing 
the plane in nearly two halves (Fig. 4b (ii)), which is an 
expected trend due to the interaction between the low energy 
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fluid collected near the inner surface (convex surface) and 
self-generated radial pressure gradient due to centrifugal force. 
In other words, these vortical motions may be the result of the 
imbalance of pressure driven secondary motions which forces 
the low-energy fluid from concave surface to convex surface 
along the circular wall. To maintain the continuity, the flow 
close to the convex surface moves towards the concave 
surface. Consequently, contra rotating vortical motions are 
developed near top and bottom surface. The vortices occupy 
almost l/3rd portion of the exit cross-section (Fig. 4b (iii)). It 
is also seen that at the exit the vortices in the vortical pair are 
not equal in strength and size. There is slight reduction in the 
cross-flow velocity vectors at the exit in comparison to the 
plane just downstream of the inflexion. 

Turbulence intensity distribution for uniform flow 

Normalized turbulence intensity distribution across 
different sections of the duct is presented in Fig. 4c.From Fig. 
4c (i), it can be seen that the turbulence intensity at inlet is low 
and fairly uniform. The turbulence intensity continuously 
increases at successive cross-sections along the duct length. 
As the flow moves into the diffuser the distortion of 
turbulence intensity profiles is minimized close to outer 
surface (Fig. 4c (ii)) of the first bend due to the bulk 
movement of fluid towards the outer surface under the effect 
of centrifugal force and self imposed pressure gradient. The 
intensity continuously increases further (Fig. 4c (iii)) close to 
outer surface of the second bend. The turbulence intensity at 
the exit of the duct is quite low (- 22.069j/kg) for major 
portion of the cross -section. The high turbulence intensity can 
be seen at the exit to be concentrated at the central portion of 
the duct. 




^22.5^22.5*1 
Fig. 5a. Normalized longitudinal velocity distribution with the 
uniform flow 
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Fig 5b. Normalized cross-flow velocity distribution with the 
uniform flow 




iU) 22.5= 22.5- 



Fig 5c.Turbulence intensity distribution with the uniform flow 
Longitudinal velocity distribution for clockwise swirl 

Figure 5 a presents the normalized longitudinal 
velocity distribution at sections 0°, 22.50/4.5° and 22.5°/22.5° 
for clockwise swirl flow at the inlet. 

At inlet, the major mass of fluid moves to the 
peripheral surface and there is fall in the center due to the 
influence of the swirler. The combined effect of centrifugal 
force due to imposed swirl flow and geometry curvature 
makes the fluid move outwards. This can also be seen from 
Fig. 5a (ii). As the flow moves downstream, the same sense of 
rotation is observed until the exit (Fig. 5b (iii)). In comparison 
to the no swirl case, the bulk of flow has shifted from the 
convex wall to the top of the duct and it is comparatively more 
uniform. 

Cross-flow velocity distribution for clockwise swirl flow 



Normalized cross-flow velocity distribution at the 
same sections as above for clockwise swirl flow at inlet is 
shown in Fig. 5b. From Fig. 5b (i), it can be seen that the 
cross-flow at inlet is completely rotational and no vortex pair 
is formed as seen for uniform flow at inlet. When the flow 
moves into the duct, there is a slight increase in the magnitude 
of cross-flow velocity in the top portion of the duct, whereas 
in the bottom portion there is a reduction in cross -flow 
velocity. This can be observed from Fig. 5b (ii). In addition, it 
is also noticed at the inlet that the cross -flow velocity core is 
slightly shifted outwards due to the downstream curvature. 
The values of the cross-flow velocities decrease in successive 
downstream sections. The center of the swirling core slightly 
shifts towards the concave surface in the first bend. The same 
sense of rotation continues in the second bend. The effect of 
swirl continues to be observed until the diffuser exit (Fig. 5a 
(iii)). In the second bend, due to a change in direction of 
curvature, the cross-flow fluid core starts moving towards 
concave surface of the second bend (Fig. 5a (iii)). From the 
cross-flow velocity distribution, it is seen that the intensity of 
inlet swirl is high enough not to allow the formation of the 
pair of vortices as seen for no swirl case. The cross-flow 
velocity shows clear circular motion across the planes in the 
clockwise direct 

Turbulence intensity distribution for clockwise swirling 

The corresponding normalized turbulence intensity 
variation at the same sections for clockwise swirl flow at inlet 
is shown in Fig. 5c. At inlet (Fig. 5c (i)), intensity values are 
higher than those observed for uniform flow velocity (Fig. 4c 
(i). However, as the flow moves downstream, the intensity 
magnitudes undergo a gradual reduction, and finally at the exit 
(Fig. 5c (iii)) the values are much lower than those observed 
for uniform flow velocity (Fig. 4c (iii)). This may be due to 
suppression of the secondary flows as a result of interaction 
between swirl flow and centrifugal force. 




Fig 5a.Normalized longitudinal velocity distribution with 
clockwise swirl 

Fig 5b. Normalized cross-flow velocity distribution with 
clockwise swirl 
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Fig 5b. Normalized cross-flow velocity distribution with 
clockwise swirl 




Fig 5c.Turbulence intensity distribution with clockwise swirl 
flow 

Longitudinal Velocity Distribution for Anti- Clockwise 
Swirl 

The normalized longitudinal velocity distributions for 
anti-clockwise swirl flow at inlet are shown in Fig 6a (i). The 
following observations and conclusions can be drawn from 
the respective figures: 

a) On comparison with clockwise swirl flow there is no 

major qualitative difference in velocity distribution. 

The minor deviations observed is the bulk flow at the 
exit for anti-clockwise has moved towards the bottom instead 
of top and the velocity magnitude is maximum for anti- 
clockwise as compared to clockwise swirl flow. This implies 
that clockwise swirl flow at inlet leads to a more uniform flow 
at exit. 

Cross-flow velocity distribution for Anti -Clockwise Swirl 

The cross-flow distributions for anti-clockwise swirl 
flow at inlet are shown in 6i. The following observations and 
conclusions can be drawn from the respective figures: 
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a) On comparison with clockwise swirl flow there is no 
major qualitative difference in cross-flow velocity 
distribution. 

The minor deviations observed is that cross-flow velocity core 
is slightly shifted towards the concave wall where as for 
clockwise it was towards convex wall. 

Cross-flow velocity distribution for Anti-Clockwise Swirl 

The cross-flow distributions for anti -clockwise swirl flow at 
inlet are shown in 6i. The following observations and 
conclusions can be drawn from the respective figures: 

a) On comparison with clockwise swirl flow there is no 

major qualitative difference in cross-flow velocity 

distribution. 

The minor deviations observed is that cross-flow velocity core 
is slightly shifted towards the concave wall where as for 
clockwise it was towards convex wall. 

Turbulence Intensity Distribution for Anti-Clockwise 
Swirl 

The turbulence intensity distributions for anti- 
clockwise swirl flow at inlet are shown in Fig. 6c. The 
following observations and conclusions can be drawn from 
the respective figures: 

a) On comparison with clockwise swirl flow there is no 

major qualitative difference in turbulence intensity 

distribution. 

The minor deviations observed is that turbulence intensity 
at the exit is relatively higher for anti-clockwise swirl 
compared to clockwise swirl but is still lower than that for 
uniform flow. 




Fig. 6a. Normalized longitudinal velocity distribution with 
anti-clockwise swirl 
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Fig 6b. Normalized cross-flow velocity distribution with anti- 
clockwise swirl 




— — >»;::.: : ::5 : 

Fig 6c.Turbulence intensity distribution with anti- 
clockwise swirl 
CONCLUSION 



On the basis of analysis and results presented, the following 
conclusion can be drawn: 

• The overall static pressure recovery increases for 
swirl flow at inlet irrespective of direction of swirl 
(clockwise or anticlockwise), the increase being 
around 40%. The loss coefficient also increases for 
swirl flow. 

• The flow distribution at the exit plane is more 
uniform for clockwise swirl in comparison to anti- 
clockwise swirl and uniform flow at inlet. 

• The secondary flow is present throughout the 
diffuser for uniform flow. A pair of vortices is 
observed at all the sections with the intensity of the 
two vortices being not equal. For swirl flow, only 
circular motion is present and no pair of vortices is 
formed. 



• Turbulence intensity variation indicates the presence 
of high turbulence intensity at the planes just before 
and after the inflexion plane. Turbulence intensity 
variation at the exit for the clockwise swirl flow is 
much lower compared to the uniform flow at the 
inlet. 
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